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Understanding the effects of geometry and rotation on 
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ABSTRACT 

Wc have developed a method to compute the possible distribution of radio emission 
regions in a typical pulsar magnetosphcre, taking into account the viewing geometry 
and rotational effects of the neutron star. Our method can estimate the emission 
altitude and the radius of curvature of particle trajectory as a function of rotation 
phase for a given inclination angle, impact angle, spin-period, Lorentz factor, field line 
constant and the observation frequency. Further, using curvature radiation as the basic 
emission mechanism, wc simulate the radio intensity profiles that would be observed 
from a given distribution of emission regions, for different values of radio frequency 
and Lorentz factor. We show clearly that rotation effects can introduce significant 
asymmetries into the observed radio profiles. We investigate the dependency of profile 
features on various pulsar parameters. We find that the radiation from a given ring of 
field lines can be seen over a large range of pulse longitudes, originating at different 
altitudes, with varying spectral intensity. Preferred heights of emission along discrete 
sets of field lines are required to reproduce realistic pulsar profiles, and we illustrate 
this for a known pulsar. Finally, we show how our model provides feasible explanations 
for the origin of core emission, and also for one-sided cones which have been observed 
in some pulsars. 
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1 INTRODUCTION 

Of the various aspects relevant for solving the unresolved 
problem of radio emission from pulsars, two are probably 
the most significant : the actual mechanism of the emission 
itself, which is still not fully understood (e.g. Zhang 2006, 
Melrose 2006) ; and the effects of viewing geometry and pul- 
sar rotation, which can significantly alter the properties of 
the pulsar profiles that the observer finally samples. The 
latter aspect has received significant attention in the recent 
years, but much still remains to be investigated. Blaskeiwicz 
et al. (1991, hereafter BCW91) were the first to work out the 
basic effects of rotation, and showed that the observed asym- 
metry between the leading and trailing parts of pulsar radio 
profiles can be due to rotation effects. Further improvements 
were carried out by Hibschman & Arons (2001), who ana- 
lyzed the first order effects of rotation on the polarization 
angle sweep. Later, Peyman & Gangadhara (2002), adapting 
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the method of BCW91, refined the formulation and showed 
that the asymmetries due to rotation can be ascribed to the 
differences in the radius of curvature of the particle trajecto- 
ries on the leading and trailing sides of the magnetic axis. By 
analysing the pulse profiles of some selected pulsars, which 
clearly show the core-cone structure in the emission beam, 
Gangadhara & Gupta (2001, hereafter GG01) and Gupta & 
Gangadhara (2003, hereafter GG03) showed that the asym- 
metry in the locations of the conal components around the 
central core component can be interpreted in terms of aber- 
ration and retardation (A/R) effects (combined effects of 
rotation and geometry) , leading to useful estimates of emis- 
sion heights of the conal components. Further refinements of 
these concepts have been carried out by Dyks et al. (2004), 
and Gangadhara (2004 & 2005, hereafter G04 & G05), Dyks 
(2008) and Dyks et al. (2009). 

All of the above said works have established that ro- 
tation effects are of significant importance in understand- 
ing the observed emission profiles of radio pulsars. What 
has been found wanting is a detailed, quantitative treat- 
ment that couples the rotation effects in the pulsar magne- 
tosphere to the possible emission physics and to the emission 
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and viewing geometries, to produce observable radio profiles. 
Some impediments to this have been recently overcome by 
Thomas and Gangadhara (2007, hereafter TG07), who have 
considered in detail the dynamics of relativistic charged par- 
ticles in the radio emission region, and obtained analytical 
expressions for the particle trajectory and it's radius of cur- 
vature. 

In this paper, we describe a scheme for simulation of 
pulsar profiles that encompasses a detailed treatment of all 
the effects mentioned above. We start with describing the 
background and motivation for the work (*J2]), then go on 
to the profile simulation method (Sj3j. We describe the main 
results from our study and their dependence on pulsar pa- 
rameters in |4] and discuss how realistic pulsar profiles may 
be obtained from our model. We also address the issues of 
core emission, one-sided or partial cones, and extension of 
our method to other models of emission physics. Our final 
conclusions are summarized in SJS] 



2 SIGNIFICANCE OF GEOMETRY AND 
ROTATION EFFECTS 

The charged particles produced in the pulsar magnetosphere 
are initially accelerated in a region very close to the polar 
cap, due to the electric fields generated by the rotating mag- 
netic field (Sturrock f971; Ruderman & Sutherland 1975; 
Harding & Muslimov 1998). After crossing the initial accel- 
eration region they enter the radio-emission domain where 
the parallel component of the electric field is screened by the 
pair plasma, and henceforth they move 'force-free'. In the ro- 
tating frame, the charged particles are constrained to move 
along the field lines of the super-strong magnetic field, the 
geometry of which is believed to be predominantly dipolar 
in the radio emission region (e.g., Xilouris et al. 1996; Ki- 
jack & Gil 1997) - the multi-polar components of the pulsar 
magnetic field are expected to be limited to much lower al- 
titudes close to the stellar surface. The accelerated charges 
are believed to produce coherent radio emission at specific 
altitudes in the magnetosphere, by a mechanism that is as 
yet not fully understood (e.g. Ginzburg et al. 1969; Ginzburg 
& Zheleznyakov 1975; Melrose 1992a, 1992b; Melrose 2006). 

Though the trajectory of the particles in the co-rotating 
frame is identical to that of the field line they are associated 
with, the trajectory will be significantly different in the ob- 
server's frame, due to the effect of rotation (TG07). The 
velocity vector of the particles will be offset from the tan- 
gent vector of the field line on which they are constrained 
to move. The value of the azimuthal angle of this offset is 
termed as the aberration phase shift <5</> a ber which depends 
on the emission altitude and also on inclination angle a and 
impact angle j3 (G05). One consequence of this rotation ef- 
fect is that the radius of curvature p of the particle trajectory 
becomes significantly different from that of the field line, as 
shown clearly by TG07. This can be intuitively understood 
as follows : on the leading side, the induced curvature due to 
rotation has the same sense as the curvature of the field lines 
and hence the net curvature of the particle trajectory in the 
observers frame gets enhanced, resulting in a reduced value 
of p, as compared to that for the corresponding field line. 
On the trailing side, the curvature of the field lines and the 
induced curvature due to rotation are in opposite directions, 



and hence they counter-act to result in a reduced effective 
curvature (or a larger value of p) for the particle trajectory. 
It was shown in TG07 that the disparity between the values 
of radii of curvature of field lines and that of the correspond- 
ing particle trajectory could be substantial. For example, for 
a pulsar with a = 90° and spin period P = 1 sec, at an emis- 
sion altitude of 0.04 of the light cylinder radius ()"l), along 
a field line with field line constant r c = 50, the ratio of the p 
estimated with and without rotation effects is more than 4. 
At an altitude 0.08 7"l, the ratio is more than 5 (see Fig. 6 in 
TG07). These ratios steeply increase for inner field lines, and 
indicate that the rotation induces substantial and unavoid- 
able differences on curvature radii and intensity of emission 
between leading and trailing sides. Further, it was shown 
in TG07 that the maximum value of the p for the particle 
trajectory, corresponding to field lines either very close to 
magnetic axis or the ones falling in the meridional plane, 
including the effects of rotation can be p max = rL/(2sina). 
While the maximum of p for these field lines in the non- 
rotating case can attain any value up to infinity. Thus, a 
proper estimation of p for the particle trajectory needs to 
take into account the effects of rotation. Unfortunately, this 
consideration has been missing in several works where p has 
been presumed to be identical to that of the dipolar field 
lines (e.g. Cheng & Zhang 1996; Lyutikov et al 1999; Gil 
et al. 2004). 

The effects of rotation and pulsar geometry can combine 
to produce interesting results in observed pulsar profiles, the 
details of which can depend somewhat on the specific mod- 
els for the emission process, including that for the distribu- 
tion of regions producing accelerated charged particles on 
the polar cap. For our immediate purposes, we use the ba- 
sic model of nested cones of emission, along with a possible 
central core, as postulated and demonstrated by several au- 
thors (e.g., Rankin 1983, 1993a, 1993b; Mitra & Deshpande 
1999; GG01; GG03). In such models, the source of emission 
could be in the from concentric rings of sparks produced 
in the polar vacuum gap, and circulating around the mag- 
netic axis (e.g., Gil & Krawczyk 1997, Deshpande & Rankin 
1999, Gil & Sendyk 2003). Each ring or cone is represented 
by a narrow annulus of field lines characterised by a defi- 
nite value of the field line constant r c appearing in the field 
line equation r = r c sin 2 9, where r is the radial distance to 
an arbitrary point on the field line and 6 is the magnetic 
co-latitude. The set of field lines can also be identified by 
5l, the distance from the magnetic axis, of the point where 
the field line pierces the neutron star surface, normalized to 
that for the last open field line (GG01). 

In the above model, if we assume that the charged par- 
ticles along a given conal set of field lines emit radiation at 
a given height in the co-rotating magnetosphere, then the 
simplest effect of rotation and pulsar geometry is to shift 
this radiation by <50 a ber in pulse longitude on the leading 
and trailing sides of the profile (with respect to the mag- 
netic meridian). This is because rotation causes the emis- 
sion beam to be offset from the local field line tangent in 
the observer's frame, hence causing the corresponding emis- 
sion component to be advanced in azimuthal phase by the 
same amount. This is the well understood aberration ef- 
fect which, in combination with the retardation effect, leads 
to phase asymmetry between the leading and trailing side 
components associated with a given cone of emission, and 



has been explored in detail by several authors (e.g., GG01; 
GG03; Dyks & Harding 2004; Dyks et al. 2004; G05). 

Furthermore, if we assume that the emission mechanism 
is such that the curvature of the particle trajectory plays 
an important role in the generation of the radio waves, as 
would be the case for models related to the curvature radia- 
tion mechanism, then rotation effects can produce significant 
changes in the strength of the intensity profiles on the lead- 
ing and trailing sides. This has been indicated by several 
authors (e.g., BCW91; Peyman & Gangadhara 2002; TG07; 
Dyks 2009). For example, TG07 have shown quantitatively 
that the ratio of the total intensity estimated for the same 
field line with and without rotation effects at an emission al- 
titude of 0.04 tl is more than a factor of 13, and at 0.08 tl it 
is more than a factor of 30 (see Fig. 8 in TG07). In extreme 
cases, this effect could lead to almost one-sided intensity 
profiles. That such effects are seen in profiles of known pul- 
sars (e.g., Lyne & Manchester 1988, hereafter LM88) is a 
strong indicator of the importance of rotation effects. 

All this motivates the importance for a method that 
can estimate the properties of the received emission, after 
including the effects described above. The necessity of such a 
3D method is pointed out by Wang et al. (2006) . A good way 
to proceed is to simulate the emission properties for a given 
choice of pulsar parameters and predict the observed profiles 
that would be seen for a range of the parameter space. These 
can then be compared with realistic pulsar profiles to gain a 
better understanding of the physical processes involved. We 
have developed such a scheme, which is described in detail 
in the following sections. 



3 PROFILE SIMULATION STUDIES 
3.1 Basic concepts 

The ultra-relativistic particles that are constrained to move 
along the co-rotating magnetic field lines suffer acceleration 
and hence emit beamed radiation within a narrow angular 
width of 2/7, centered on the direction of the instantaneous 
velocity vector, v. This emission is aligned with the local 
tangent vector b to the field line in the co-rotating reference 
frame. However, as mentioned earlier, in the observer's ref- 
erence frame, the velocity vector is offset from the field line 
tangent vector (GG01, G05). Furthermore, the radius of cur- 
vature of the particle trajectory will be different from that of 
the associated field lines (TG07). To begin with, for a given 
pulsar geometry we choose a ring of field lines specified by 
a single value of field line constant r c (or by the equivalent 
value of Sl) and look for all possible emission spots along the 
field line that can contribute in the observer's line-of-sight 
direction. For the observer to receive significant radiation, 
we impose the condition that the unit velocity vector, v, 
should align with the unit vector along the line-of-sight, n, 
such that 

n - v = 1 . (1) 

For any given pulse phase of observation, we find the pos- 
sible emission spots on the specified ring of field lines that 
meet this criterion. For each emission spot, we compute the 
emission altitude, r, the emission angles 9 and 4> (defined 
in sec. i)3.2[) . and the radius of curvature of the particle 
trajectory, p. 
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Figure 1. A schematic picture that depicts the geometry asso- 
ciated with the emission region. The XYZ frame is co-rotating 
with the magnetic axis rh around spin axis f2, and the X'Y'Z' 
frame is the non-rotating frame of the observer. 

We then couple the basic curvature radiation model to 
the above picture, using the value of p in order to obtain 
estimates of the specific intensity that would be seen by 
the observer at a given pulse longitude. Of the several emis- 
sion mechanisms proposed for pulsar emission, the curvature 
emission model is perhaps the most natural and favoured one 
(Gil et al. 2004). However, as we argue later, our method 
is flexible enough to incorporate other viable variants of 
pulsar emission models. Using the curvature radiation for- 
mulation, we compute the intensity that would be seen at 
any observing frequency, u), and not just for the charac- 
teristic frequency, u c . We sweep the line of sight through 
discrete rotation phases and, for any given field line, search 
for the points that contribute at a given to, and calculate 
the emission parameters like emission height, radius of cur- 
vature and, finally, the intensity of the radio emission re- 
ceived. This method allows us to generate a super-set of all 
possible emission profiles that would be observed at a given 
frequency. 

3.2 Details of the method 

Any emission point in the pulsar magnetosphere can be lo- 
cated by the coordinates r, 9 and in a coordinate system 
where the Z-axis is parallel to the magnetic axis m, and the 
XZ-plane is the plane containing the magnetic and rotation 
axes (see Fig. [I]). This coordinate system can be called as 
the magnetic coordinates, where 9 and <fi are the magnetic 
co-latitude and the magnetic azimuth, respectively. This co- 
ordinate system co-rotates with the pulsar. Another coor- 
dinate system, identified by X'Y'Z' (the observer's frame), 
can be defined such that the line-of-sight vector is parallel 
to the X' Z' -plane containing the magnetic axis and the ro- 
tation axis at rotation phase cj>' = 0, and designated as the 
meridional plane M. The Z'-axis is parallel to the rotation 
axis f2. The X'Z'— plane makes an azimuthal angle <f>' with 
M during rotation. 

In the X'Y'Z' coordinate system, the location of the 
emission spot can be specified by the values of r, 9' and 
<f)p (see qXp . The expressions for 9 and <j) are given in G04, 
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and the expression for 9' as a function of the pulsar rotation 
phase cf>' is given in G05. A point of emission on a dipolar 
field line can be expressed as r = r c sin 2 9. When rotation 
effects are not included, it turns out to be a trivial exercise 
to trace emission spots based on the expressions for r(t), 
9 and (j>. For a given a, /3 and Sl combination, the value 
of r(<j>') — r c sin 2 [6 (((>')], the radius of curvature p can be 
readily found (Eq. 4 in G04). However, the estimation of r 
and p corresponding to the emission spot becomes difficult 
when the effects of rotation are included. 

As mentioned before, when rotation is invoked, the ob- 
server will receive peak radiation when the particle velocity 
vector v (rather than the tangent to the field line) becomes 
parallel to the line-of-sight h, where n = [sin (, 0, cos Q, and 
f = a + p. The total velocity of the charged particle will be 
the vector sum of the component parallel to the magnetic 
field, and the component in the azimuthal direction due to 
the co-rotation of the field (TG07). The component in the 
azimuthal direction will make the total velocity offset from 
the direction of the local field tangent. The analytical solu- 
tion for the radial position of the trajectory of the charged 
particle, including the rotation effects, is derived in TG07. 
We employ the zeroth order solution that gives : 



r{t) = — — cn(A - f2 m0 1) 

S'mO 



(2) 



where f2 m o = sin a, and cn(z) is the Jacobi Cosine function 
(Abramowitz & Stegun f972), A is a constant, t is the affine 
time and c is the speed of light. 

For a given a, /3 and a rotation phase (/>' , the radial 
distance for an emission spot can, in principle, be found by 
solving Eq. JT}. The rotation phase </>' is defined as the pro- 
jected angle on the azimuthal plane between n and m (see 
Eq. f| A 1 1 1> in §A). The above-said analytical expressions for 
the coordinates of the accelerated charged particle are to be 
as invoked as functions of affine time. The aberration phase 
shift <50 a ber needs to be known in advance for solving Eq. {TJ 
and finding the emission altitude, since the angles 9 and <j> 
are aberrated due to the rotation (see G05 for details, in- 
cluding the analytical expression for 5</> a ber)- Nonetheless, 
the emission altitude is apriori required for estimating the 
#</>aber, thus making the problem non- linear. An analytical 
solution for Eq. {TJ is nearly impossible owing to the bulky 
trigonometric terms. Approximations like <5<^ a b cr ~ r/rt,, 
will severely limit the expected precision of the solution. A 
straight forward numerical solution for Eq. JTJ) also encoun- 
ters more or less similar difficulties owing to the aforesaid 
reasons. Hence we have developed special algorithms which 
are applicable for solving Eq. |T]) under such conditions. We 
have devised an 'exact' method and also an 'approximate' 
method to compute the location of possible emission spots, 
which are suited for specific parameter regimes. These are 
described below. 



bo. But when the effects of rotation are invoked, the emission 
beam at Po gets mis-aligned with bo and goes out of the line- 
of-sight. Hence the radiation from Po will not be received by 
the observer. However, another emission spot Pi on the same 
ring of field lines can have the v parallel to n and contribute 
emission in the direction of the observer. Hence the observer 
will receive radiation from Pi provided Eq. ([T]) is satisfied 
for Pi. Let the aberration phase shift at Pi be <5</> a ber; then 
the updated values 9 — 9(<f)' + 5^> a ber) and <j> = <?H0' + 50 a ber) 
at the point Pi so that the emission is aligned with the line 
of sight to the observer (G05). 

Hence the basic idea that is invoked in the computation 
of the exact method can be briefed as: for a given <f>' , begin 
with the value of the emission height, which is estimated 
in the non-rotating case, and find the trial value of 50 a ber- 
Then solve Eq. ([TJ numerically to find an improved value of 
emission height. Continue the iteration till the solution for 
Eq. |TJ satisfactorily converges. The main steps of algorithm 
are briefly described below : 

(i) Choose a specific combination of a, f3, 7, r e and to, 
and a fixed rotation phase cj>' . 

(ii) Make the first estimate for the aberration angle <50 a b O r 
using the trial input values of 9(<f)'), </>(</>') and r, which fol- 
lows from r = r c sin 2 9((f>'). 

(iii) Using the 54>ahei estimated above, the angles <j> and 
9 are re-calculated with the rotation phase incremented 
by <50 a bcr- Henceforth update: (f>((f>') — > (f>(<j>' + <50 a ber) and 
9(<t>')->9(4>' + 5^ hcl ). 



(iv) Estimate v and hence the unit vector v ; 
the angles <j> and 9, found in step (iii) 



v/|v| with 



(v) Estimate the affine time 't' that satisfies the matching 
condition n • v = 1. Hence find the improved value of r(t). 

(vi) Recalculate <5</> a bor with r = r(t) and repeat the cal- 



culation from step (iii) till convergence is achieved for 't'. 

(vii) Using the improved value of r(t) find v, a, and p. 

(viii) Using the p estimate the spectral intensity. 

(ix) Find the angle ?7 m i s = cos - (n • v). 

We choose to call this method as 'exact method' since 
the computation employs the exact expressions for the rel- 
evant quantities. The explicit expressions for magnetic co- 
latitude 9 and magnetic azimuth (f> are given by Eq. (|A19l l 
and Eq. (| A20|l . respectively, in §A. The expressions for veloc- 
ity v and acceleration a are given by Eq. (|A1() and Eq. (|A2|I 
respectively, and the expression for radius of curvature p is 
given in Eq. (| A22|l . Sample results from this method are 
shown in the figures in SjC] The angle r; m i s gives the residual 
difference between the line-of-sight and the estimated v at 
the end of the iterations. It's ideal value is zero and hence, 
the final residual value obtained is a measure of the preci- 
sion of the solution : a smaller value of ?7 m i s indicates a more 
precise determination of the emission spot. 



3.2.1 The 'exact' method 

For a ring of field lines specified by a field line constant r e 
and for a given rotation phase </>', we consider a point Po on 
a field line such that the unit-vector of the local field line 
tangent bo is parallel to n in the non-rotating case. In the 
absence of rotation, the emission beam from the accelerated 
particles moving along the field line should be aligned with 



3.2.2 Alternative or 'approximate' method 

We have devised an alternative or approximate method, for 
the estimation of r, p and the related quantities, in cases 
where we encounter 'extreme' values of parameters. Such 
regimes are often combination of large values of a, very low 
values of /3, and field lines close to magnetic axis (Sl < 0.5). 
The exact method encounters difficulties for such regimes in 
that the numerical solutions of Eq. ([T]) for affine time 't' 



Rotation effects on pulsar intensity profiles 5 



often do not give satisfactory convergence. So we resort to 
an approximate method that is suitable for this regime. By 
this method, we expect to determine the emission height 
and the radius of curvature with comparable precision to 
the exact one, for leading and the trailing parts of the pulsar 
profile. The estimates of this method have been optimised by 
comparison with the estimates of the aforesaid exact method 
in a common parameter regime where both the methods give 
reliable results. The scheme of the approximate method is: 
for a given rotation phase <j>' , first calculate r, 9 and cj> in 
the non-rotatmg case, and then use it to find approximate 
values of r, 6 and cfi in the rotating case. The details are 
provided in ijB] 



3.3 Computing the intensity profiles 

As explained earlier, we sweep the line of sight through dis- 
crete rotation phases, and using the afore-said steps given 
in £|3 . 2 . II and i]3.2.2l we find the parameters like emission 
altitude r, and the radius of curvature of the particle trajec- 
tory, p. Then we estimate the spectral intensity for a given 
frequency, uj, for particles for a given Lorentz factor, 7, by 
using the standard curvature radiation formula (e.g. Jackson 
1972) : 

dl/duj = V3 — 7 u I K 5 / 3 (u')du , (3) 

c Ju 

where u = uj/uj c and the characteristic frequency uj c = 
1.5 7 3 e/p. According to Eq. (J3}, the spectral-intensity curve 
should peak at u = 0.286, which corresponds to p/p P = 1, 
where the parameter p p can be defined as 

p p = 0.286 x I.57 3 - . (4) 

uj 

Invoking this parameter helps in easy identification of the 
peak points in the spectral-intensity plots. 

An important feature of this method is that it computes 
the contribution to the observed intensity for any frequency 
uj, different from uj c . Often in literature, significant contri- 
bution of intensity to the observer is presumed to be con- 
centrated near the characteristic frequency uj c = 1.57 3 c/p, 
thus providing an in- built frequency selection criteria (e.g., 
Melrose 2006). However, the spectral intensity curve for cur- 
vature radiation has non-negligible amount of power emitted 
at a significant range of frequencies different from uj c . Our 
formulation thus allows for a more complete treatment of 
the amount of emitted intensity and its reception by the 
observer. 

Since our present analysis necessitates only the com- 
putation of relative intensities of the simulated profiles, in- 
voking a single particle emission model for the curvature 
emission do not alter the results in a significant manner. 
The high luminosity of the pulsars demands imposing co- 
herence on the emission, perhaps in the form of bunched 
emitting sources, and the process behind the formation of 
such bunches is still being investigated. In a simple manner, 
coherent emission from a bunch with charge Q can be al- 
ternatively expressed as the emission from a single particle 
with the same charge Q. So, the relative intensities are not 
affected by this simplification and hence considering single 
particle emission do not tamper with the physics behind the 



emission. Further discussions regarding this factor will be 
followed in later sections. 

3.4 Typical outputs 

We have computed the parameters of emission in the mag- 
netosphere by implementing the method described above. 
The free parameters are the following : a, f3 (pulsar geome- 
try), fi (pulsar rotation frequency), 7 (Lorentz factor of the 
particles), Sl (field line location) and uj (radio frequency 
of observations). For the sake of brevity of presentation, we 
give results for a single fixed value of SI = 2tv (i.e. a spin 
period of 1 sec), and for a relatively narrow range of rota- 
tion phases of about —10° to +10° around zero (fiducial) 
phase. As is discussed later, frequency turns out to be a rel- 
atively weak parameter in comparison to other strong ones 
that influence profile evolution. Hence, for the simplicity of 
analysis, we have restricted the frequency to a single value of 
610 MHz. For a chosen pulsar geometry, the emission loca- 
tions are estimated for each discrete rotation phase and for 
a set of discrete choices of Sl. Further, the specific intensity 
values are estimated for a set of discrete values of 7, for a 
fixed value of uj. 

The typical outputs are shown in the figures in iJC] 
Fig. IC II shows the basic outputs from a typical simulation 
run for estimating the location of emission regions, for a 
fixed pulsar geometry, for a set of Sl values (0.1, 0.3, 0.5 
and 0.7). The following quantities are shown, as function 
of rotation phase, in separate panels for each Sl value : 
the estimated emission altitude, r; the computed radius of 
curvature of the particle trajectory, p; the <f> and 6 for the 
emission spots in the magnetosphere; and the mis-alignment 
angle, r\mis, which is a measure of the accuracy of the re- 
sults. Fig. IC2l shows the computed intensity profiles for each 
choice of Sl values (in separate panels) for a set of 7 val- 
ues (200, 300, 400, 600, 1000, 1500), for a fixed observation 
frequency of 610 MHz. These two figures illustrate the ba- 
sic results. The effect of varying pulsar geometry can now 
be explored to understand the variety of intensity profiles 
that are possible. Most important outputs are r/ri,, p/r^ 
and specific intensity plots, and the successive figures show 
only these quantities. 

Figures [C3I to [C7l then investigate the effects as a func- 
tion of pulsar geometry (i.e. different combinations of a and 
P values), with the range of Sl values and 7 values confined 
to those illustrated in Figs. IC 1 1 and IC2I 



4 RESULTS AND DISCUSSION 
4.1 General results 

We first discuss the general results and trends that are de- 
duced from our simulation studies, as illustrated in the re- 
sults displayed in the plots given in ijC] 

4-1.1 Emission heights 

The heights for the allowed emission spots have a minimum 
value near the magnetic meridian (</>' = 0), with smoothly 
increasing values on the leading and trailing sides. However, 
the variation of height with rotation phase is asymmetric 
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such that the increment with rotation phase is always faster 
for the trailing side than for the leading side. Whereas this 
increase with rotation phase is purely geometric, the asym- 
metry in this is due to the modification of particle trajec- 
tories produced by rotation. On the leading side, the emis- 
sion beam bends in the direction favourable to rotation and 
hence advances in azimuthal phase by 54>'a heT from that of 
the corresponding field line tangent. Hence at a fixed <j>' , 
an emission spot located at a lower emission height than in 
the non-rotating case will satisfy Eq[T] In contrast, on the 
trailing side, the bending of the emission beam causes a lag 
in azimuthal phase by <5^4ber from the corresponding field 
line tangent. This lag can be compensated by an emission 
spot located at a different altitude than the non-rotating 
case. Hence an emission spot located at a different emis- 
sion height will satisfy Eq[T]and contribute radiation to the 
observer. 

Both the value of the minimum height (at </>' = 0), as 
well as the asymmetry, are larger for the inner field lines as 
compared to the outer field lines. This supports the intuitive 
expectation that outer field lines would be visible out to 
much larger pulse longitude ranges than inner field lines. 
Further, it is seen that the minimum height and asymmetry 
of the variation increase with geometry, being more for larger 
values of a and j3. Also, pulsars with larger values of a will 
have a larger range of variation of allowed emission heights, 
for the same field line. It is reasonable to surmise that a 
similar value of emission height can be seen recursively for 
several combinations of larger and smaller values of a and 
Sl, for a certain cj>'. Thus we find that the values of a and 
Sl dominantly decide the range of emission heights. Even 
otherwise, the a dependence of the emission height can be 
directly understood from Eq.©, since the expression r(t) is 
explicitly dependent on the value of ct. 



4-1.2 Radius of curvature 

The radius of curvature inferred for the possible emission 
spots also varies signficantly with rotation phase, and it is 
significantly asymmetric between leading and trailing sides. 
If rotation effects were not considered, this radius of curva- 
ture would be same as that of the corresponding field line 
and be symmetric around zero pulse phase. Given that the 
observed radius of curvature of the particle trajectory is a 
combination of the curvature of the field line and curvature 
introduced due to rotation, the observed asymmetry is a 
rotation effect and can be understood as a combination of 
two effects: first, as the allowed heights are different on the 
leading and trailing sides (for the same phase ±<j>' on either 
side of the zero phase), the radius of curvature of the field 
line itself would be different; second, on the leading side, the 
curvature of the field line gets combined with that induced 
by rotation (TG07), resulting in a lower value of the p for 
the particle trajectory. Whereas on the trailing side, the two 
curvatures are opposed and hence the net curvature is re- 
duced, leading to larger values of p. In fact, for some field 
lines, rotation induced curvature can cancel the curvature 
due to the field lines, at some points on the trailing side, 
resulting in sharply peaked curves for p for certain values of 
5l, as seen in the plots in Figs. lC4llC5l fc [C6l The variation 
of p on the leading side is more or less steady, while on the 



trailing side it often varies very rapidly due to the aforesaid 
reasons. 

The trend is that the asymmetry seen in p will be higher 
for larger values of a and smaller values of Sl and hence 
it is a combined effect of both. As mentioned earlier, the 
rotation effects are more at higher a and hence the larger 
asymmetry. Since the range of emission altitudes covered by 
the emission spots for inner field lines are higher than that 
of the outer field lines, the corresponding values of p also 
will have a higher range and a higher asymmetry than the 
ones for the outer field lines. However, there are variations 
in the amount of asymmetry with in the field lines when the 
a varies, and hence a steady variation of asymmetry with a 
will not be observed as for the emission heights. 

Comparing the p/ph plots for figs. IC5I fc IC61 we can 
find on the leading side that the radius of curvature gets 
reduced when a increases from 60° to 90° for all field lines. 
But on the trailing side the behaviour is slightly different. 
The inner field lines (Sl = 0.1 & 0.3) have the radius of 
curvature slightly reduced on the trailing side, while an in- 
crement in radius of curvature is seen for outer field lines 
(Sl = 0.1 & 0.3), when a increases from 60° to 90°. An- 
other comparison of the p/ph plots for the figs. IC1I & IC5I 
obviously shows the same trend for the leading side. How- 
ever, the variation of p/ph on trailing side shows a slightly 
different behaviour, varying among field lines. The variation 
of p on the trailing side is not in a steady pattern as on 
the leading side. The competing curvatures due to rotation 
induced curvature and the intrinsic field line curvature gives 
a highly varying pattern for p on the trailing side. 

4-1-3 Spectral Intensity, 1^ 

The derived spectral intensity curves reflect the asymmetry 
inherited from the variation of p with pulse phase, combined 
with the effects of the value of 7. In particular, it is read- 
ily seen that the intensity dramatically evolves with 7. For 
lower values of 7, the I u has a stronger leading part while 
for higher values of 7 the J„ has a stronger trailing part. 
This effect can be better understood by considering the pa- 
rameter p p (defined in Eq. ©), which gives the value of 
p at which the spectral intensity peaks, for given values of 
frequency and 7. For values of p greater than or less than 
p p , the spectral intensity falls monotonically. Further, the 
peak value of the spectral intensity also depends on the spe- 
cific value of gamma, as per Eq. ©. Hence, the variation of 
spectral intensity with pulse longitude can be inferred from 
that of p with longitude, for different values of gamma. This 
is illustrated in fig. IC2I where p/p p and the corresponding 
Iuj^ 1 ) are plotted side by side as a function of the rotation 
phase, for specific combinations of parameters. 

Three different cases are useful to consider. For situ- 
ations where p/p P is greater than 1.0 for the entire pulse 
window, the spectral intensity curve shows a maximum at 
<j> = and falls asymmetrically on either side, with the re- 
duction in intensity being larger for the trailing side, due to 
the faster increase of p. This effect, which is seen for rela- 
tively small values of 7 (less than 400-600), naturally leads to 
asymmetric pulse profiles, with possibilities for sharply one- 
sided profiles. It is interesting to note that in some cases, 
the intensity on the trailing side can drop to negligible val- 
ues, compared to its value at the corresponding longitude 
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on the leading side. This could be a natural explanation for 
the one-sided cones reported in literature, and is discussed 
in more detail in i )4.6l 

For situations where p/p P is less than 1.0 for the entire 
pulse window, the spectral intensity curve shows a minimum 
at <f>' = and rises asymmetrically on either side, with the 
increase being larger on the trailing side. However, the con- 
trast in the intensity levels between leading and trailing sides 
is typically not as high as for the first kind. This behaviour 
is seen for relatively large values of 7. 

For intermediate cases, where values of p/p P less and 
greater than 1.0 can occur at different points in the pulse 
longitude window, we see more complicated shapes for the 
spectral intensity curves, including multiple maxima at dif- 
ferent pulse longitudes. 

For a given geometry and sl value, the transition 
through these 3 different cases can take place as 7 is var- 
ied over a range of values. Thus, lower values of 7 tend to 
produce profiles with strong leading and weak trailing in- 
tensities which get converted to weak leading and strong 
trailing kind as the 7 increases to very high values (e.g. 2nd 
and 3rd panels from the top in fig. IC2|) . Though the con- 
trast in intensity is less for the latter, the absolute value of 
the spectral intensity is higher, due to the 7 dependence in 
Eq. ©. 

Most of the asymmetric intensity profile effects become 
more dramatic for inner field lines and for more orthogonal 
rotators (larger a) and smaller values of j3. 

4-1-4 9 an d 4> 

The values of 6 and <j) are asymmetric on the leading and the 
trailing sides of the profiles, while they are symmetric in the 
non-rotating case (G04). This asymmetry is also an effect of 
rotation. Since 6 and (j> are functions of <f>' their values are 
affected by the aberration phase shift which is different on 
leading and the trailing sides. Their values are dominantly 
decided by a and /3. As expected, the shape of the <f> curve 
closely resembles the S-shape of the polarization angle curve 
(see fig. [Cl]|. 

4-1.5 Mis-alignment angle, ?7 m is 

The Mis-alignment angle ry m i s , defined as ?7 m is = cos _1 (n-v), 
gives an estimate of the offset between n and the estimated 
v. In principle, for a perfect estimation of the emission point, 
the line of sight should exactly coincide with the velocity 
vector, and hence ry m i s should be zero. However, in actual 
computations, 7) m i s always has a small, finite value. A quick 
check of the accuracy of the computation is provided by the 
value of ?7 m i s : a lower value implies a higher precision of 
estimation of the emission spot, and vice versa for a higher 
value. A rough classification that can be taken for the pre- 
cision of the estimation is: a value of r/mia <C 1 indicates 
a highly precise estimation of the emission spot, and vice 
versa for ?7 m i s 3> 1. By this scheme, we find that there is 
satisfactory precision for all estimations for r within 20 % of 
tl (see fig. IC8[) . In some cases the r; m i s exceeds 1, but only 
when r/rL > 0.2. However, according to established obser- 
vational results radio emission heights are restricted within 
10 % of tl for normal pulsars (e.g. Kijak 2001), and hence 
our method is quite satisfactory in this regime of interest. 



4.2 Effects of Parameters 

The above described behaviour of the height of emission 
spots, radius of curvature and spectral intensity are strongly 
dependent on the parameters like geometry, field line loca- 
tion, radio frequency and Lorentz factor of the particles. In 
some cases, there is a complex interplay between the depen- 
dencies on these different parameters. Here, we explore some 
of these effects in detail. 

The generic effects of 7, a and /3, u) and Sl are listed 
briefly below in an order that may characterize the hierarchy 
of their effects on total intensity profiles. 



4-2.1 Inclination angle, a 

The parameter a is a major driver of the effects of rotation, 
and has the strongest influence on our results and conclu- 
sions. The rotation effects (leading-trailing asymmetry of r, 
p and Iu,) are more prominent for large values of a and less 
for small values of a. The range of emission altitudes is found 
to be relatively high for lower values of a, and relatively low 
for higher values of a, being the lowest for a = 90°. Like 
wise, p appears to reach higher values for higher a and vice 
versa for lower a. 



4-2.2 Normalized foot value of the field lines, Sl 

The effect of moving from inner to outer regions of the mag- 
netosphere (increasing Sl values) also has a very dramatic 
effect on the results. Rotation effects are strongest for the 
innermost field lines, and decrease significantly for larger Sl 
values. For relatively small values of Sl (usually < 0.3), the 
leading part has emission heights that vary relatively gently 
with </>', whereas the trailing part shows steeply rising emis- 
sion height. The emission heights become less asymmetrical 
for increasing values of Sl- The values of p steadily increase 
with decreasing Sl (i.e. inner field lines) on the trailing side. 
Dramatic effects such as very large, peaked values of p on 
the trailing side, owing to the mutual cancellation of intrin- 
sic and rotation induced curvatures, are seen only on inner 
field lines. This peak shifts closer towards zero pulse phase 
as the value of Sl becomes smaller. For outer field lines, the 
profiles are much more symmetric, and since p/p p < 1 for 
a significant range of pulse phase on either side of 4>' = 0, 
the profiles more often exhibit minima at (j)' = 0, even for 
moderate values of 7. 



4-2.3 The lorentz factor, 7 

The effect of 7 has a very clear signature on the asymmetry 
of the spectral intensity profiles. For lower values of 7, there 
can be strong asymmetries with leading side stronger than 
the trailing side, and maxima at <f)' = 0. For larger values 
of 7, the sense of this asymmetry can reverse, with trailing 
side becoming stronger than the leading, and a minima at 
4> = ; however, the degree of the asymmetry, as measured 
by the ratio of the intensities at corresponding longitudes, 
is generally less than that for the case for the low 7 values. 
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4-2.4 Emission frequency, to 

The frequency of emission, uj, acts as a counter to the ef- 
fect of 7, though in a relatively weak manner, as can be 
understood from Eqs. @ and Q. Thus, an increase in 
oj produces changes which can be compensated by a corre- 
sponding change in 7 by a factor proportional to to 1 / 3 . In 
certain cases, this could result in profiles where the sense of 
asymmetry between leading and trailing sides could reverse 
over a large enough range of radio frequencies. Such effects 
are seen sometimes in some real profiles. 



4.3 Realistic profiles 

One of the significant results from our simulation studies 
is that the possible regions of emission associated with a 
given annular ring of field lines (characterised by a constant 
value of Sl) are visible over a wide range of pulse phase, 
albeit with different intensity levels. This aspect, combined 
with the results for field lines with different values of Sl, 
leads to the conclusion that a very large fraction of the pul- 
sar magnetosphere is potentially visible to us. This results 
in simulated pulsar profiles that are very different from the 
observed profiles of real pulsars which appear to have well 
defined emission components, restricted in pulse phase ex- 
tent to occupy only some fraction of the on-pulse window. 
This disparity with the observed profiles persists even after 
we incorporate into our simulations the models of discrete, 
annular conal rings of emission. 

Hence, in order to reproduce realistic profiles matching 
with observations, we need some additional constraints for 
the emission regions. In the most general case, such non- 
uniformities in the distribution of emission regions can exist 
in any of the three coordinate directions, viz. r, 9 and cj>. 
Non-uniformity of emission in the 6 direction is achieved in 
some sense at the basic level, by considering only discrete 
sets of Sl values for active emission regions. As discussed, 
this is not enough to constrain the intensity variations to 
reproduce realistic pulsar profiles. 

The possibility of non-uniform emission along the (f> co- 
ordinate could help produce discrete emission components 
in the observed profiles. As seen in fig. ICI I for a given value 
of Sl, the emission at different pulse longitudes originates at 
widely different <j) locations. If the sources of charged parti- 
cles were located only at fixed (f> points along the ring of con- 
stant Sl, then these could be arranged to modulate the sim- 
ulated intensity pattern with a suitable "window" function, 
to obtain discrete emission components in the observed pro- 
file. However, the well known phenomenon of sub-pulse drift 
argues against this being a viable option. Sub-pulse drift, 
which is now believed to be fairly common in known pul- 
sars (e.g. Velterwede et al. 2006), wipes out any azimuthal 
discretization of sources of emission in the pulsar magneto- 
sphere - it would lead to "filling up" of the intensity average 
profile over a given range of pulse longitude, as is observed 
in drifting sub-pulses that occur under any discrete emission 
component in known pulsars. 

The third option is to have non-uniform emission in 
the radial direction, i.e. preferred heights of emission for a 
given set of field lines. Since the contributions at different 
pulse phases are from different heights, this would naturally 
lead to non-uniform distribution of intensity in pulse phase, 
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Figure 2. PSR B2111+46 at 333 MHz : The emission height for 
each Sl value associated with a particular component is plotted 
in panel (a). The simulated and un-modulated spectral intensity 
curve corresponding to each Sl value is plotted in panel (b). The 
simulated sub-components, after applying Eq. (0 for the best fit 
values given in Tablc[l]arc shown in panel (c). The sum-total of 
the simulated sub-components giving the final profile (solid curve) 
is shown superposed with the observed profile (dotted curve) in 
panel (d). 



resulting in realistic looking pulse profiles. The idea of pre- 
ferred heights of emission in the pulsar magnetosphere is not 
entirely new - the "radius to frequency mapping" model for 
pulsar emission postulates different heights for different fre- 
quencies, with the height of emission increasing for decreas- 
ing frequency values (Kijak & Gil 1997,1998; GG01, GG03). 

Preferred emission heights of emission with a spread 
in the r direction, can be modeled as a multiplication of 
the spectral intensity I u (</>') with a modulating function, 
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Table 1. The parameters for simulating the profiles of PSR B2111+46 4 - 4 Profiles for PSR B2111+46 : a test case 



Frequency 


i" 


Hi 


AHl 


K 


Sl 


Lorentz factor 


MHz 




Km 


Km 






7 


333 


1 


1500 


600 


1.8 


0.08 


750 




2 


1834 


500 


0.14 


0.18 


750 




3 


3800 


500 


0.7 


0.3 


500 


408 


1 


300 


850 


2.2 


0.09 


750 




2 


1200 


750 


0.1 


0.22 


750 




3 


3000 


400 


0.9 


0.35 


500 


610 


1 


200 


700 


1.7 


0.13 


750 




2 


500 


450 


0.3 


0.31 


700 




3 


2600 


500 


1.5 


0.35 


550 



1 i = 1 represents the core component, 

i = 2 represents the inner conal component, 

I = 3 represents the outer conal component. 



Fu{4>')- For modelling a profile as a sum total of emissions 
from a core and several discrete conal regions, the modulated 
spectral intensity can be expressed as 







Ft, 



Al exp 



'■<>( 



2AHL 



(5) 
(6) 



where the index i represents, a corresponding pair of leading- 
trailing components presumed to be arising from a particu- 
lar ring of field lines; while i = 1 exclusively represents the 
central core component. Here N is the total number of such 
discrete emission components (for example, N = 3 would 
correspond to a 5 component profile forming a central core 
and two pairs of conal components), I u {(j>') is the total spec- 
tral intensity from all field lines combined, while II, ((/)') is 
the spectral intensity from the ith ring of field lines. For a 
given emission region, HI, represents the mean height, and 
AHl, represents the spread of the region. The variable rj) (</>') 
represents the values of emission altitude for each value of 
</>' estimated by the simulation method as described earlier, 
for the ith ring of field lines. 

To map this intensity as a function of rotation phase 
as seen in the observer's frame, the effects of retardation 
and aberration need to be included explicitly. The effect of 
aberration is estimated by default and the <j>' is inclusive 
of the aberration phase shift. The retardation phase shift 
is to be estimated from the value of r corresponding to the 
emission spot. The rotation phase corresponding to I u is 
updated after adding the retardation phase shift dcj>' Tet with 
<j}' , and this is represented by the mapping of the ordered 
pair (cf)',Iu,) -» (4>' + d(f)' T et>Itu)- The d<j>' Tet can be estimated 
as (G05) : 

#ret = —(r-n), 

where r = re r and the expression for e r is given in [JB] The 
height of emission (HI,) and the normalized foot value (S L ) 
corresponds to the peak of the ith component of the profile. 
The AHl, and A* u are model parameters. 



Using the afore said methods for simulation of pulse pro- 
files, we have attempted to reproduce the intensity profiles 
of PSR B2111+46 obtained from EPN data base and GMRT 
data, at multiple radio frequencies. This pulsar has a multi- 
component profile, with a well identified core component 
and 2 cones of emission (e.g. Zhang et al. 2007). ft has a 
rotation period of 1.014 sec and a = 14 and /3 — —1.4 (Mi- 
tra & Li 2004). The other parameters used in the simulation 
are listed in Table \T\ The values of emission heights HI, and 
field line locations S L for the discrete emission components 
are the values from estimates employing the method given 
in Thomas and Gangadhara (2009). The zero phase of the 
profile is fixed on the basis of the analysis of core emission 
of this pulsar, using the method developed in the the same 
work. 

The simulation method is illustrated in fig. [2] Values of 
Sl corresponding to the core and conal components are em- 
ployed in generating the emission height plots in panel (a) of 
fig. [2] The corresponding spectral intensity plot for each of 
these Sl values, for the final best fit choice of 7 (in Table[T]) 
is shown in panel (b) of this figure. The individual compo- 
nents generated after applying the best fit height function 
are shown in panel (c) and the sum total intensity profile 
is shown in panel (d), along with the observed profile. Best 
fits of these profiles to the observed data were obtained by 
varying Al, and AHl, in the function F£, (</>'), and by varying 
the value of 7 in the range 100 to 1000. The same procedure 
is repeated for 408 MHz and 610 MHz profiles and the re- 
sults are shown in fig. [3] All the final parameters and best-fit 
results are summarised in Table [T] 

An encouraging first order match between the simu- 
lated and observed profiles has been achieved (see panel (d) 
of fig. [2] & panels (b) and (d) of fig. [3}. The core compo- 
nent is quite well fit for most of the cases, and so are the 
leading conal components. There is some mismatch in the 
widths of the conal components, especially for the trailing 
side, where the real data shows a smoother blending of the 
components, compared to the simulated profile where the 
components appear more narrow and relatively well sepa- 
rated. It is remarkable that with a single value of Al, for a 
leading-trailing pair of cones of emission, the ratios of the 
peak values of the intensity of the leading and trailing com- 
ponents of the cones match so well with the real data. It 
is also interesting to note that the best fit values for 7 are 
very similar for a given emission component, at different fre- 
quencies, supporting a model of a common bunch of accel- 
erated charged particles being responsible for the emission 
at different frequencies. Further, that the spread of 7 values 
across the different emission components is also quite small, 
indicates very similar operating conditions over most of the 
magnetosphere. The best fit values for AH^, , though reason- 
able, are somewhat large in amplitude, indicating somewhat 
extended emission regions in the magnetosphere. 

We note that these relatively large values of AH^, and 
some of the limitations of the fits may be due to the lack of 
some generalizations in our model. These include factors like 
coherency of emission, a realistic spread of 7 values around 
the mean values obtained here, as well as a realistic spread 
in the values of Sl due to finite thickness of the rings of 
emission on the polar cap. Whereas a detailed treatment of 
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all of these is beyond the scope of this work and will be 
taken up later, some basic inferences can still be drawn. For 
example, if a small range of Sl values around the mean is 
considered, it is easy to argue that much of the width of a 
profile component can be filled up by radiation from such 
a bunch of field lines. This can be understood from panel 
(a) of fig. [2] where a line of constant height intersects the 
curve for a given field line at two points, one each on the 
leading and trailing side. The phase of this point of intersec- 
tion will move systematically as we go to neighbouring field 
lines. Hence, wider profile components can be achieved with 
smaller values of AH^. Furthermore, due to the asymmetry 
in the emission height curves, the shift in phase with change 
of S L is more on the trailing side, which would naturally 
lead to broader component widths and better 'blending' of 
the components in the profile, something that is not as eas- 
ily achievable by having a large range of emission heights 
(as the shift of phase for a given separation of heights on a 
given field line is lesser on the trailing side). One indicator 
of the significance of the spread of S L values is the amount 
by which needs to be changed to move the peak of one 
conal component to the point half-way to the peak of the 
next conal component. Not very surprisingly, our rough es- 
timates show that the required change in S L is close to the 
half-way point to the S L value of the next cone, which would 
indicate a closely packed structure of concentric rings. 

The component profiles could be further influenced by 
considering a distribution of 7 values associated with the 
emitting particles. We have found that significant shifts in 
the peaks of the leading and trailing pair of components for 
simulated profiles are obtained for lower 7 values (7 < 500) , 



while the peak positions appear almost frozen for increasing 
7 values. Thus it is realistic to assume that a spread of 7 
values can broaden the emission components. This factor 
also may reduce the AHl, required to effect a good fit. 

Nevertheless, we would like to point out that there is 
only one unique combination of the parameters that can 
produce a profile which is similar to the observed one. We 
have not found any degenerate combination of values for the 
parameters that are shown in TablefT] Thus, the similarity 
of the simulated profiles with the observed ones gives an 
assurance that, we should be able to simulate the observed 
profiles with greater similitude with a model overcoming the 
above-said limitations. 



4.5 Core emission 

The generation of the profile components for PSR B2111+46 
described above naturally leads to a discussion on the core 
emission. In fact, the study of the phenomena of core emis- 
sion has spawned enormous amount of literature. Perhaps 
the most notable ones are the landmark work by Rankin 
(1983) that systematized the pulsar emission profile into 
'core' and 'cone', and the succeeding works by Rankin 
(1993a & 1993b) that further developed the core-cone clas- 
sification scheme. The hollow cone model was invoked to ex- 
plain the geometry (e.g. Taylor & Stinebring 1986) and the 
origin of core emission. Radhakrisnan and Rankin (1990) 
have conjectured that the emission mechanism for cores 
might be different from that of cones, owing to the behaviour 
of polarization position angle curve near the core being dif- 
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ferent from the rotating vector model. However, there are 
no satisfactory theoretical grounds for postulating diverse 
mechanisms for cores and cones. A major difficulty that cur- 
vature radiation encounters in explaining the core emission 
is the insufficient curvature of the almost straight field lines 
in a region relatively close to the pulsar polar cap. Since the 
intensity of emission is proportional to 1/p 2 , the values of p 
provided by the intrinsic curvature of the field lines is too 
large and hence insufficient to generate enough intensity of 
emission typically observed for core component. This fac- 
tor even prompted invoking other emission mechanisms for 
explaining core emission (e.g. Wang et al. 1989). 

In our simulation studies, the presence of the core com- 
ponent comes about quite naturally. It can be seen from all 
the plots of spectral intensity in iJC] that the emission from 
regions near the profile centre is comparable (for higher 7 
values) or even somewhat higher (for lower 7 values) than 
that from regions in the wings of the profile. This happens 
because we get low values of p for inner field lines near 
</>' ~ 0, which are comparable to that of outer field lines, 
and this occurs consistently for all the combinations of a 
and ft (see the panels for p/rx in jJCj. The reason is that 
rotation induces significant curvature into the trajectory of 
particles, even though they are confined to move along the 
nearly straight inner field lines. 

The forces of constraint act in such a way that the par- 
ticle is hardly allowed to deviate away from the field line on 
which it is moving, and the resulting scenario is discussed in 
detail in TG07. Due to the co-rotation of the field lines and 
the action of the aforesaid forces of constraint the charged 
particles are added with a velocity component in the direc- 
tion of rotation, which is nearly perpendicular to the velocity 
component parallel to the field line, in the observer's frame. 
This induces an additional curvature in the trajectory of 
the particle and makes it significantly different from that of 
the field line curvature in the observer's frame of reference 
(TG07). Hence the trajectory of charge particles moving on 
almost straight field lines near the magntic axis can have a 
highly curved trajectory and hence a relatively low value of 
radius of curvature that is significantly different from that of 
the field lines. This scenario allows for significant emission 
near the central region of the profiles. By applying Eq. 
and Eq. ([6]) appropriately, as described earlier, profile shapes 
resembling strong core components can be easily generated. 
Hence by applying our method, we provide a natural expla- 
nation for core emission, that circumvents the issue of too 
high p that precludes a strong core component with curva- 
ture emission. 

In the simulation of the profiles for PSR B2111+46, we 
find that the core originates from have relatively inner field 
lines and lower emission heights than the cones. Assuming 
the same mechanism of emission, viz. curvature radiation, 
for the core and conal component, we are able to produce a 
simulated core component that matches quite well with the 
observed profile. We notice that the best-fit values for the 
amplification factor found in the simulation for the core 
component (Table [T]) are comparable to those of the cones. 
These values are not unduly high, considering the situation 
that the density of plasma should be relatively higher for the 
lower altitude and hence an additional factor for relatively 
stronger emission at lower altitudes. Our profile-matching 



of PSR B2111+46 thus shows that strong core emission can 
originate from inner field lines due to curvature emission. 

4.6 Partial cones 

According to LM88, partial cone profiles are the ones in 
which one side of a double component conal profile is either 
missing or significantly suppressed. These are recognised by 
the characteristic that the steepest gradient of the polariza- 
tion position angle is observed towards one edge of the total 
intensity profile, instead of being located more centrally in 
the profile as the rotating vector model postulates. LM88 
speculated that this happens when the polar cap is only 
partially (and asymmetrically) active. It is significant that, 
out of the 32 pulsars listed in LM88 that display the partial 
cone phenomena, as many as 22 have the steepest gradient 
point occurring in the trailing part of the profile. In other 
words, most of the partial cone profiles show a strong leading 
component and an almost absent trailing component. 

There are two possible scenarios that have been postu- 
lated to explain partial cones : (1) only a part of the polar 
cap is active (this works for both kinds of partial cones) or 
(2) the A/R effects are so large as to shift the entire active 
region of the intensity profile towards the leading side (this 
works for the strong leading type partial cones, which are 
the majority). However, Mitra et al (2007) studied several 
pulsars with partial cones with very high sensitivity obser- 
vations and found that the almost-absent parts of the cones 
do flare up occasionally and show emission for about a few 
percentage of the total time. This tends to rule out both 
the scenarios above, and requires an explanation where the 
intensity is naturally suppressed in one side of the cone. 

In our simulation studies, one sided cones appear as a 
natural by-product. We notice that for smaller values of a, 
inner field lines and lower 7 values, the intensity profile is 
almost always significantly suppressed on the trailing side, 
as compared to the leading side. The reason for this is quite 
obvious. As explained earlier, the p for inner field lines is 
highly asymmetrical between the leading and trailing sides 
- it remains more or less steady on the leading side, while on 
the trailing side it shoots up to a high value and then falls. 
Whenever the p shoots up such that p/p p ^> 1, the spec- 
tral intensity is significantly reduced. On the other hand, 
on the leading side we mostly have p/p P ~ 1 and hence the 
spectral intensity is significant there. For relatively lower 
values of 7, p p reduces and hence there is a greater chance 
of having p/p p 3> 1, while for higher values of 7, p/p P drops 
down and eventually becomes closer to 1. Hence the inten- 
sity plots shown in [JC] are stronger on the leading side at 
low 7, and stronger on the trailing side at high 7. How- 
ever, it is to be noted that (i) the values of 7 required to 
achieve stronger trailing side profiles are very high - usu- 
ally significantly more than 1000; whereas, for more typical 
values of 7, we get the stronger leading side profiles and 
(ii) the intensity contrast obtained for the stronger leading 
side profiles is much larger and striking, compared to that 
for the stronger trailing side profiles. Both these facts argue 
naturally for a strong preponderance of one sided cones with 
stronger leading side profiles, as is statistically seen in the 
results of LM88. 

To further illustrate the idea, we have generated profiles 
as shown in fig. IC7I that resemble partial cone profiles by 
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Table 2. The parameter values employed in fig |C7l 



Panel a S L HI AH^ 

No. [°] [°] Km Km 

1 90 1 0.3 1000 200 

2 90 1 0.5 1200 150 

3 60 1 0.3 1000 200 

4 60 1 0.5 1400 150 

5 30 1 0.1 2500 500 

6 30 1 0.3 1000 150 



using our simulation technique for specific combinations of 
parameters, which are listed in Table [2] The thin line curve 
represents the un-modulated profile, which is simulated by 
assuming that the emission is uniform all along the field 
line, while the thick line shows the final modulated profile. 
The active region of the final profile is clearly shifted to the 
leading side, due to the afore said behaviour of p. The sup- 
pression of the intensity on the trailing part in comparison 
to the leading side is seen in all the plots and is most dra- 
matic for the inner field lines, for low values of 7, and for 
large values of a. 

4.7 Studying the mechanisms of emission : future 
prospects 

In this section, discuss some of the future possibilities from 
the present work. Though we have developed a model un- 
der certain specific conditions and demonstrated some useful 
results from the same, it has significant potential for appli- 
cability under diverse circumstances and conditions. The r 
and p of the emission spot are two of the fundamental in- 
gredients for computing the intensity of emission within any 
model of radio emission for pulsars. These values, along with 
other parameters that we have calculated in our method af- 
ter explicitly taking into account of effects of rotation and 
geometry, are applicable for any model of radiation that pre- 
cepts the condition embodied in Eq.(JT}, i.e. having the ra- 
diation beam aligned with the velocity vector and line-of- 
sight. Thus the present method of computing r and p is well 
suited for studying curvature radiation models in vacuum 
approximation. The profiles simulated by these models can 
be compared with the observed ones to check their veracity. 

Though we have employed single particle curvature ra- 
diation formulation, it is well known that this cannot explain 
the extremely high luminosities seen in typical pulsar radio 
emission. Coherent emission from bunches of charged par- 
ticles have been argued to be necessary for explaining the 
high luminosities (eg. Ginzburg et al. 1969, Melrose 92, Mel- 
rose 2006). The model of coherent emission constructed by 
Buschauer and Benford (1976) considered relativistic charge 
and current perturbations propagating through the bunches 
with N number of charges, which boosted the emission much 
above typical N 2 factor. Further, they have shown that the 
characteristic frequency will be significantly shifted to higher 



values than the typical ~ 1.57 3 c/p. However, the extremely 
short lifetime of these moving sheets of plasma (bunches) 
made it implausible to radiate, and due to this reason these 
emission models were almost forgotten. In later years, the 
possibility of formation of Langmuir micro-structures (soli- 
tons) due to the collective behavior of the plasma brought 
back the possibility of bunched radiation (Asseo 1993). It 
was shown that the radiation from such a bunch could be 
expressed by just using the classical formula for curvature 
radiation (Asseo 1993). Melikidze et al. (2000) considered 
the three component structure of charge distribution for soli- 
tons in the pulsar magnetosphere and obtained a different 
spectral intensity distribution from that of the classical for- 
mula for curvature radiation. However, Gil et al. (2004) used 
single charged bunches of charge Q as equivalent to a single 
particle with the same charge Q, to explore the effects of the 
surrounding plasma on the curvature emission, and showed 
that sufficient luminosity could be produced from curvature 
emission that matches with the observed luminosity of pul- 
sars. 

As mentioned in the earlier H3.3I our estimates and re- 
sults corresponding to altitude, radius of curvature, mag- 
netic azimuth and magnetic colatitude are equally valid for 
the case of coherent and incoherent emission, as long as 
the vacuum approximation is invoked. This is because the 
peak of the emitted beam will be aligned with the direc- 
tion of velocity for emission from a source moving at ultra- 
relativistic speeds, de-facto in vacuum approximation. Hence 
the premise contained in Eq.fl} for the computation of these 
quantities will remain valid for both of the cases. For the case 
of a simple model of coherence for a bunch of net charge Q, 
the spectral intensity profile estimated will be similar to that 
of the emission from a single particle with charge Q, and 
likewise the relative intensity will also be the same. Hence 
the results that we have drawn upon spectral intensity are 
valid for the simple case of coherent emission too. However, 
invoking models of coherent radiation with additional fea- 
tures apart from a simple coherent model, may push the 
intensity estimates to significantly different values and the 
resulting shape of the intensity profile will be considerably 
altered. Two such examples are mentioned in the following. 

Buschauer and Benford (1976) has shown that both 
intensity profile and characteristic frequency will be al- 
tered if the allowance is made for the propagation of 
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a charge and current density wave through the coherent 
bunch. Considering this model we find that it can alter the 
shape of the computed spectral intensity curve correspond- 
ing to a given field line, from that of the present results. 
This is mainly because of the reason that characteristic fre- 
quency u c will be shifted to a higher value than in the case 
of single particle curvature radiation. Another case is the 
spectral intensity formula for emission from solitons having 
a three component charge structure (Eq.(12) in Melikhidze 
et al. 2000) which also will yield significantly different esti- 
amtes for spectral intensity, from that of the spectral inten- 
sity estimated for the single particle emission. Both of these 
models are treated in the vacuum approximation and hence 
they satisfy the condition embodied in Eq.JT]), i.e. having the 
radiation beam aligned with the velocity vector and line-of- 
sight. This ensures that the method of estimation and hence 
the results corresponding to altitude, radius of curvature, 
magnetic azimuth, magnetic colatitude etc. will be applica- 
ble for these two cases also. The only quantity that is altered 
by the inclusion of these models, from a single particle case, 
is the spectral intensity estimate. Nevertheless, these models 
can be quite easily incorporated into our simulation studies, 
simply by modifying the form of the spectral intensity ex- 
pression that is used. 

In the emission models where the effects of the sur- 
rounding plasma are considered, the peak of the radiation 
beam may be offset from the velocity vector by a finite 
angle (Gil et al. 2004). This requires a modification to the 
condition in Eq.fTJ) such that n • v = J7 max where ?7 m ax is 
the value of the angle of offset by which the peak of the 
emission beam is offset from the velocity vector. Coupling 
this with some modifications to our method can deliver 
the values of r and p appropriate for this case too. The 
analysis and results that ensue from all of the above said 
considerations will be discussed in our forthcoming works. 



5 SUMMARY 

We have developed a method to compute the probable lo- 
cations of emission regions in a pulsar magnetosphere that 
will be visible at different pulse longitudes of the observed 
profile. The effects of geometry and rotation of the pulsar 
are accounted in a detailed manner in this method, which is 
a very useful new development. Our method includes 'exact 
' and 'approximate' techniques for carrying out the estima- 
tion of the relevant emission parameters. The 'approximate' 
method is useful for certain extreme regimes of parameter 
space, and for faster computation of the results. The mis- 
alignment angle, which provides a good check of the accu- 
racy of the computations, shows that our method achieves 
satisfactory precision. Besides the exact location of possible 
emission regions, we are able to compute several other use- 
ful parameters like the height of emission, and the radius of 
the particle trajectory at the emission spot, the azimuthal 
location of the associated field line etc., for different com- 
binations of pulsar parameters like a and /3. Further, using 
the classical curvature radiation as the basic emission mech- 
anism (which is apt for a debut level analysis), we are able 
to compute the spectral intensity from any emission spot. 
By assuming a uniform emission all along the field lines, we 



have estimated the spectral intensity for a range of pulse 
phase that the line of sight sweeps through. We have dis- 
cussed how realistic looking pulsar profiles can be generated 
from these generalized intensity curves, by assuming spe- 
cific range of emission heights along specific rings of field 
lines. We have illustrated the capabilities of these methods 
by generating simulated profiles for the test case of the pul- 
sar PSR B2111+46, and have shown that fairly good match 
with observed profiles can be achieved. We have also shown 
how further detailed (and practical) considerations can help 
improve this match. We have shown how our results offer a 
direct and natural explanation for the puzzling phenomena 
of partial cones that are seen in some pulsar profiles. Our 
simulations also provide a direct insight into the generation 
of the core component of pulsar beams. Finally, we have 
indicated how our method can be extended to incorporate 
more sophisticated models for the emission mechanism and 
produce intensity profiles for the same. These, as well as ex- 
tension to polarized intensity profiles, will be taken up as 
future extensions of the work reported here. 
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APPENDIX A: VELOCITY, ACCELERATION 
AND RADIUS OF CURVATURE OF PARTICLE 
TRAJECTORY 

Al Expression for velocity and acceleration 

The expressions for velocity v, acceleration a, radius of cur- 
vature p, etc. that are used in the computation described 
in 33T2T1 and 3X231 are provided here (see §C in TG07 for 
details). The velocity and acceleration (in spherical polar 
coordinates) of the charged particle in the laboratory frame 
can be denned as . 

dr „ d9' 

V = ^ er + r ^T 

and 



eg 



rsme— e,, 



(Al) 



d v 
~~dt 
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The corresponding unit-vectors and their derivatives are 
given by 

sin 9' (cos (f)' p x + sin <f)' p y ) + cos 9' z , (A3) 
cos #'(cos 4>p x + sin <f>' p y) — sin#'z , (A4) 
— sin 4>' p x + cos cj)' p y , (A5) 
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(A6) 
(A7) 
(A8) 



Here x, y and z denotes the unit vectors along the X' , Y' and 
Z' axes as described in £13.21 Using the relation r = r e sin 2 9 
valid for a point on a static dipolar field line, the following 
derivatives are found out: 

d9 d9 dr tan 9 dr , , . „ N 

(A9) 



2r dt 



and 



(A10) 



dt dr dt 

d£ d0_dB_ 
~dt ~ ~~d~9 dt 

The angle (j>' p is the azimuthal phase between the radial vec- 
tor to the position of the charged particle, and the fidu- 
cial plane containing the line-of-sight and the rotation axis. 



difference between the line-of-sight and the magnetic axis 
and it can be defined as 



(f)' — cos -1 (fix • fhj_), 

nj_ = n — z(n ■ z), 

m± — rh — z(rhz), 

rh = {sin a cos <f>' , sin a sin <f> , cos a}, 

n = {sin (, 0, cos £} 



(All) 
(A12) 
(A13) 
(A14) 
(A15) 



A<j) is the azimuthal phase difference between the the ra- 
dial vector to the position of the charged particle, and the 
magnetic axis and it is given as 



A(f) — cos 1 ^ 



cos 9 sin a + cos a cos <f> sin f 



(A16) 



A2 Expressions for 9, <f> and <50 a b C r 

The following expressions which are used in the computation 
are given in G04 and G05. 



r = cos 1 [cos a cos C, + sin a sin ( cos 0'] 



= - cos 
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= tan 
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- | cos 
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lV8 + cos 2 T - sin 2 r 
sin £ sin <f>' 



(A17) 
, (A18) 



cos £ sin a — cos a sin £ cos cf>' 
— cos - 1 [cos a cos 9 — sin a sin 9 cos 4>] 
The aberration phase shift 50 a b ol - is given as (G05) 
r sin 9' cos G 



(A19) 
(A20) 



tan £ cot ip + 



(A21) 



tl sin £ sin ip 
where the angles Q and ip defined in G05. 

A3 Expressions for p and b 

The radius of curvature is found out using the expression 
(see TG07 for details) 



P = T^- r (A22) 
v x a| 

The expressions for v and a are given in Eq. (|A1[1 and 
Eq. (|A2|l . The position vector of an arbitrary point on a field 
line in the coordinate system-XyZ, with the Z-axis point- 
ing in the direction of rh is given by 



r = r c {sin 3 9 cos <f>, sin 3 9 sin cj>, sin 2 9 cos 9}, (A23) 



A- r 



(A24) 



Hence 



b' ± A<f>. The angle <j)' is the azimuthal phase 



and A = R ■ I is the product of / (Inclination) and R (ro- 
tation) matrices (G04). Then the field line tangent in the 
coordinate system- A"!" Z' is given by b = dr t /d9 and 
b = b/|b|. 



APPENDIX B: DETAILS OF THE 
APPROXIMATE METHOD 

In the 'approximate method', we utilize the parameters cor- 
respoding to emission spot Po in the non-rotating case as 
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input values for estimating the emission spot Pi correspond- 
ing to the rotating-case. To facilitate this, we employ a few 
approximations to estimate the 8 and <f> corresponding to 
Pi . For a ring of field lines specified by a field line constant 
r c and for a given rotation phase <f>' , we take a point Po 
on a field line such that the unit-vector of the local field 
line tangent bo is parallel to n in the non-rotating case. 
If the effects of rotation are negliglected then the emission 
beam from the accelerated particle moving along the field 
line should be aligned with bo. When the effects of rotation 
are invoked the emission beam at Po will be aligned with 
v instead of bo which makes the bo to be offset by an az- 
imuthal angle S<j>^ (aberration phase shift at Po) with n. 
Hence the radiation from Po will not be recieved by the ob- 
server. However, another emission spot Pi on the same ring 
of field lines at a different rotation phase <f>i can have the 
v parallel to n and contribute emission in the direction of 
the observer. Let the unit-vector of the local tangent be bi 
at Pi. The observer will recieve radiation from Pi provided 
the azimuth angle between bi and bo will be equal to the 
aberration phase shift <5<^ber a ^ ^ • This inevitably leads to 
the condition for the reception of radiation from Pi as 



cos" 1 (Wo - bxi) = <5^ 1 b ) er , (BI) 
where 

b ± = b-Cl(Cl- b), (B2) 

bx = fa (B3) 

Wo = b ± (0(<l>'), </>(</>'), </>'), (B4) 

bj.i = bx(0i(0i), </>i(0'i), <&), (B5) 

0i Oi) = 0(^ + 5^), (B6) 

Mti) = <K<&+8</>+) (B7) 



and bxo and bxi are the unit vectors of projections of the 
bo and bi on the equatorial plane, respectively. The (j>i is 
the azimuthal phase between rh and n corresponding to Pi. 
The phase shifts 8(f>g and 84>'$ are neccessarily introduced for 
shifting the emission spot from Po to Pi. The exact values 
of 8cj)' g and Scf)'^ can be found out by concomitantly solving 
Eq. ()B1|I and Eq. ([T]) for the point Pi . Since the calculations 
that ensue can become very cumbersome we evade it and 
instead resort to seperate approximations appropriate for 
the leading and trailing sides. 



B2 Trailing side 

For the trailing side, the phase is assigned as (/>[ — <f>' fol- 
lowed by the approximation 

6<f>4 = -6<j>' e = S^ el . 

Here, the point Pi is found out while keeping the azimuth 
angle between rh and n unchanged, which is different from 
the method availed for the leading side. The emission alti- 
tude at Pi is readily found out as r\ = r c sin 2 [0i (<?!>' — S<f>^} )] 
at the rotation phase cj>'. 

Henceforth for the point Pi, the v and a are found 
from the values of r, 0, <f> that are estimated by the meth- 
ods given above; further the radius of curvature and spectal 
intensity are computed. It needs to be verified that the few 
approximations invoked above are not at the cost of the 
accuracy of the estimation of emission spot. This can be 
verified by calculating the angle (f] m i s ) at Pi. As mentioned 
before, the value of r) m i s should be ideally zero for a perfect 
estimation of the emission spot. A non-zero value of 7? m i s 
angle indicates a less than perfect estimation of the emis- 
sion spot. The approximate method gives reasonably pre- 
cise results (?7 m is <S 1°) with in r/rL < 0.2, but gives large 
errors (?7mis 2> 1°) if the estimated emission heights exceed 
this limit. Since the observational results confirm that radio 
emission heights for normal pulsars are limited with in 10 
% of rx, and our region of interest is restricted with in this 
range of emission heights, the precision of this method is 
satisfactory for our needs. 



BI Leading side 

For the leading side we assign cj>'i = cf>' — 5<j>^ and make 
the approximation that 

Thus we find that 0i(0i) = ((/>') and = (j>((j)'). This 

implies that the values of radial distance ro = r c sin 2 [0(0')] 
at Po and ri = r c s'm 2 [0i(4>' 1 )] at Pi are equal. Hence the 
emission spot at Pi on the leading side can be readily ob- 
tained by re-assigning the rotation phase for the ordered 
pair at P, -> Pi as (0', r ) -> ((/>[, r ) = (0' - 5(p^ er , n). 
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Figure CI. Simulation results for a = 30° and f3 = 1°, 
as a function of pulse longitude, <f>' . The first row shows 
the estimated emission height, r, in Km, and as a frac- 
tion of rjj. The second row shows the estimated radius of 
curvature, p, as a fraction of n,, and the azimuthal an- 
gle, cj>. The last row shows 6 and the mis-alignment angle 
*?mis- The results are plotted for 4 values of Sl : Sl = 
0.1 (tiny dash), 0.3 (small dash), 0.5 (medium dash), 0.7 (large dash). 



APPENDIX C: SAMPLE RESULTS 

Sample results from the simulation studies are illustrated 
with a series of figures, which are explained in detail in the 
main text. 
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Figure C2. Simulation results for a = 30° and = 1°. The 
ratio p/pp is plotted in the panels in the first column, while the 
corresponding spectral intensity (in arbitrary units) is plotted in 
the panels in the second column, for the different choices of Sl- 
Each panel has results for different choices of 7, ranging from 200 
to 1500, in varying step sizes. 
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Figure C3. Simulation results for a = 30° and = 2°. The emis- Figure C4. Simulation results for a = 30o and = 3°. Sec 

sion altitude r/r L and the radius of curvature p/r L , the spectral caption for Fig.[C3]for details, 
intensity (in arbitrary units) are plotted with respect to <f>' , for 
the same range of Sl and 7 values, as in Fig. IC2I 
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Figure C6. Simulation results for a = 90° and/3 = 1°. See the 
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Figure C7. Illustrating partial cones : The simulated intensity 
profiles (thin lines) for different 7 values and the corresponding 
modulated profiles (thick lines) are plotted. The peak of the mod- 
ulated profile touches the corresponding simulated profile. The 
rclavcnt parameters used for each panel are given in Table [2] 
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Figure C8. The angle »7 m i s for the approximate method is 
plotted with respect to <j>' , for 3 different geometries in the 3 
panesl, for different values of Sl = 0.05 (tiny dash), 0.075 (small 
dash), 0.1 (medium dash), 0.125 (large dash). Next, for column 
2, S L = 0.2 (tiny dash), 0.225 (small dash), 0.3 (medium dash), 
0.325 (large dash). 
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